Chapter 25 



Deterministic diffusion 



This is a bizzare and discordant situation. 



M.V. Berry 



(R. Artuso and P. Cvitanovic) 




he advances in the theory of dynamical systems have brought a new life to 



Boltzmann's mechanical formulation of statistical mechanics. Sinai, Ruelle 



and Bowen (SRB) have generalized Boltzmann's notion of ergodicity for a 
constant energy surface for a Hamiltonian system in equilibrium to dissipative sys- 
tems in nonequilibrium stationary states. In this more general setting the attractor 
plays the role of a constant energy surface, and the SRB measure of sect. 16.1 is 
a generalization of the Liouville measure. Such measures are purely microscopic 
and indifferent to whether the system is at equilibrium, close to equilibrium or far 
from it. "Far for equilibrium" in this context refers to systems with large devia- 
tions from Maxwell's equilibrium velocity distribution. Furthermore, the theory 
of dynamical systems has yielded new sets of microscopic dynamics formulas for 
macroscopic observables such as diffusion constants and the pressure, to which 
we turn now. 

We shall apply cycle expansions to the analysis of transport properties of 
chaotic systems. 

The resulting formulas are exact; no probabilistic assumptions are made, and 
the all correlations are taken into account by the inclusion of cycles of all periods. 
The infinite extent systems for which the periodic orbit theory yields formulas for 
diffusion and other transport coefficients are spatially periodic, the global state 
space being tiled with copies of a elementary cell. The motivation are physical 
problems such as beam defocusing in particle accelerators or chaotic behavior of 
passive tracers in 2-dimensional rotating flows, problems which can be described 
as deterministic diffusion in periodic arrays. 

In sect. 25. 1 we derive the formulas for diffusion coefficients in a simple phys- 
ical setting, the 2-dimensional periodic Lorentz gas. This system, however, is not 



509 



CHAPTER 25. DETERMINISTIC DIFFUSION 



510 



Figure 25.1: Deterministic diffusion in a finite horizon 
periodic Lorentz gas. (T. Schreiber) 



the best one to illustrate the theory, due to its complicated symbolic dynamics. 
Therefore we apply the theory first to diffusion induced by a 1 -dimensional maps 
in sect. 25.2. 



25.1 Diffusion in periodic arrays 

Chaos happens - let's make a better use of it. 
— Edward Tenner 

The 2-dimensional Lorentz gas is an infinite scatterer array in which diffusion of a 
light molecule in a gas of heavy scatterers is modeled by the motion of a point par- 
ticle in a plane bouncing off an array of reflecting disks. The Lorentz gas is called 
"gas" as one can equivalently think of it as consisting of any number of pointlike 
fast "light molecules" interacting only with the stationary "heavy molecules" and 
not among themselves. As the scatterer array is built up from only defocusing 
concave surfaces, it is a pure hyperbolic system, and one of the simplest non- 
trivial dynamical systems that exhibits deterministic diffusion, figure 25.1. We 
shall now show that the periodic Lorentz gas is amenable to a purely determin- 
istic treatment. In this class of open dynamical systems quantities characterizing 
global dynamics, such as the Lyapunov exponent, pressure and diffusion constant, 
can be computed from the dynamics restricted to the elementary cell. The method 
applies to any hyperbolic dynamical system that is aperiodic tiling M = U«er M« 
of the dynamical state space M by translates M„ of an elementary cell At, with 
T the abelian group of lattice translations. If the scattering array has further dis- 
crete symmetries, such as reflection symmetry, each elementary cell may be built 
from & fundamental domain M by the action of a discrete (not necessarily abelian) 
group G. The symbol M refers here to the full state space, i.e.,, both the spatial 
coordinates and the momenta. The spatial component of M is the complement of 
the disks in the whole space. 
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Figure 25.2: Tiling of M, a periodic lattice of reflect- 
ing disks, by the fundamental domain M. Indicated is 
an example of a global trajectory x(t) together with the 
corresponding elementary cell trajectory x(t) and the 
fundamental domain trajectory x(t). (Courtesy of J. -P. 
Eckmann) 

We shall now relate the dynamics in M to diffusive properties of the Lorentz 
gas in M. 

These concepts are best illustrated by a specific example, a Lorentz gas based 
on the hexagonal lattice Sinai billiard of figure 25.2. We distinguish two types 
of diffusive behavior; the infinite horizon case, which allows for infinite length 
flights, and the finite horizon case, where any free particle trajectory must hit a 
disk in finite time. In this chapter we shall restrict our consideration to the finite 
horizon case, with disks sufficiently large so that no infinite length free flight is 
possible. In this case the diffusion is normal, with x(t) 2 growing like t. We shall 
discuss the anomalous diffusion case in sect. 25.3. 

As we will work with three kinds of state spaces, good manners require that 
we repeat what tildes, nothings and hats atop symbols signify: 

fundamental domain, triangle in figure 25.2 

elementary cell, hexagon in figure 25.2 

full state space, lattice in figure 25.2 (25.1) 

It is convenient to define an evolution operator for each of the 3 cases of fig- 
ure 25.2. x(t) = f'(x) denotes the point in the global space M reached by the 
flow in time t. x{t) - f'(xo) denotes the corresponding flow in the elementary 
cell; the two are related by 

h t (xo) = f(x Q ) - f(xo) € T , (25.2) 

the translation of the endpoint of the global path into the elementary cell At. The 
quantity x(t) = f'(x) denotes the flow in the fundamental domain At; f'(x) is 
related to f\x) by a discrete symmetry g € G which maps x(t) e M to x(t) € M . chapter 21 

Fix a vector /3 e M. d , where d is the dimension of the state space. We will 
compute the diffusive properties of the Lorentz gas from the leading eigenvalue of 
the evolution operator (17.11) 

s{p) - \im-\og(^ m - x) ) M , (25.3) 

t—. >oo t 




( 
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where the average is over all initial points in the elementary cell, x e M. 

If all odd derivatives vanish by symmetry, there is no drift and the second 
derivatives 



JLJL 



2dD u = ——s(J3) 



= lim -((x(f) - x)i{x{t) - x)j) M 
/?=() 1 



yield a diffusion matrix. This symmetric matrix can, in general, be anisotropic 
(i.e., have d distinct eigenvalues and eigenvectors). The spatial diffusion constant 
is then given by the Einstein relation (17.13) 



1 v- d 1 

D =T^W* m 



1 , 

= ton TT<(^(0 - <?) >M 



where the i sum is restricted to the spatial components of the state space vectors 
x = (q,p), i.e., if the dynamics is Hamiltonian, the sum is over the d the degrees 
of freedom. 

We now turn to the connection between (25.3) and periodic orbits in the ele- 
mentary cell. As the full At — > At reduction is complicated by the non-abelian remark 25.5 
nature of G, we discuss only the abelian M —> M reduction. 



25.1.1 Reduction from M to M 

The key idea follows from inspection of the relation 

veAl 

| Al| = J M dx is the volume of the elementary cell At. Due to translational symme- 
try, it suffices to start with a density of trajectories defined over a single elementary 
cell At. As in sect. 17.2, we have used the identity 1 = J M dy 6(y - x(t)) to moti- 
vate the introduction of the evolution operator £'(y, x). There is a unique lattice 
translation h such that y - y — h, with the endpoint y e M translated back to the 
elementary cell, and f'(x) given by (25.2). The difference is a translation by a 
constant lattice vector n, and the Jacobian for changing integration from dy to dy 
equals unity. Therefore, and this is the main point, translation invariance can be 
used to reduce this average to the elementary cell: 

^<m-x) } _L C dxdy^ <f ' <x) - x) 6(y-f(x)). (25.4) 
\M\ J x , y eM 

As this is a translation, the Jacobian is \dy/dy\ - 1. In this way the global f\x) 
flow, infinite volume state space averages can be computed by following the flow 
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f'(xo) restricted to the compact, finite volume elementary cell M. The equation 
(25.4) suggests that we study the evolution operator 



(25.5) 



where x{t) = f'(x) e M is the displacement in the full space, but x, f'(x), y e M. 
It is straightforward to check that this operator satisfies the semigroup property 
(17.27), 



f. 

Jm 



dz£' 2 (y,z)£ n (z,x) = £ 2+n (y,x) 



For /5 — 0, the operator (25.5) is the Perron-Frobenius operator (16.10), with the 
leading eigenvalue e s ° = 1 because there is no escape from this system (see the 
flow conservation sum rule (20.17)). 

The rest is old hat. The spectrum of £ is evaluated by taking the trace 



section 18.2 



tr£ f = f dx^ x) 5{x-x{f)) . 
Jm 



Here h t (x) is the discrete lattice translation denned in (25.2). Two kinds of orbits 
periodic in the elementary cell contribute. A periodic orbit is called standing 
if it is also periodic orbit of the infinite state space dynamics, P^ p (x) = x, and it 
is called running if it corresponds to a lattice translation in the dynamics on the 
infinite state space, P^''(x) = x + h p . We recognize the shortest repeating segment 
of a running orbit as our old 'relative periodic orbit' friend from chapter 9. In the 
theory of area-preserving maps such as the standard map of example 7.7 these 
orbits are called accelerator modes, as the diffusion takes place along the momen- 
tum rather than the position coordinate. The traveled distance h p - hj {xq) is 
independent of the starting point xq, as can be easily seen by continuing the path 
periodically in M. 

The final result is the spectral determinant (19.6) 



det(s06)-J?t) - ]~[exp - 



1 e (fi-n p sT p )r 



r=l 



det(l-M£)| 



(25.6) 



or the corresponding dynamical zeta function (19.15) 



|An| 



(25.7) 
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The dynamical zeta function cycle averaging formula (20.24) for the diffusion 
constant (17.13), zero mean drift (x,) = , is given by 



D = 




1 



1 



, (-l) k+l (n pi +---+h pk ) 2 



(25.8) 



2d (T\ 



2d (T\ 



where the sum is over all distinct non-repeating combination of prime cycles. The 
derivation is standard, still the formula is strange. Diffusion is unbounded motion 
across an infinite lattice; nevertheless, the reduction to the elementary cell enables 
us to compute relevant quantities in the usual way, in terms of periodic orbits. 

A sleepy reader might protest that x{T p ) - jc(0) is manifestly equal to zero for 
a periodic orbit. That is correct; h p in the above formula refers to a displacement 
x(T p ) on the infinite periodic lattice, while p refers to closed orbit of the dynamics 
f'(x) reduced to the elementary cell, with x p a periodic point in the closed prime 
cycle p. 

Even so, this is not an obvious formula. Globally periodic orbits have x 2 = 0, 
and contribute only to the time normalization (TV. The mean square displace- 
ment (x 1 ^ gets contributions only from the periodic runaway trajectories; they 
are closed in the elementary cell, but on the periodic lattice each one grows like 
x(t) 2 = (h p /T p ) 2 t 2 = v 2 t 2 . So the orbits that contribute to the trace formulas 
and spectral determinants exhibit either ballistic transport or no transport at all: 
diffusion arises as a balance between the two kinds of motion, weighted by the 
l/|Ap| measure. If the system is not hyperbolic such weights may be abnormally 
large, with l/|A p | « 1/T p a rather than l/|A p | « e~ T > lA , where A is the Lyapunov 
exponent, and they may lead to anomalous diffusion - accelerated or slowed down 
depending on whether the probabilities of the running or the standing orbits are 
enhanced. section 25.3 

We illustrate the main idea, tracking of a globally diffusing orbit by the as- 
sociated confined orbit restricted to the elementary cell, with a class of simple 
1 -dimensional dynamical systems where all transport coefficients can be evalu- 
ated analytically. 



25.2 Diffusion induced by chains of 1 -dimensional maps 

In a typical deterministic diffusive process, trajectories originating from a given 
scatterer reach a finite set of neighboring scatterers in one bounce, and then the 
process is repeated. As was shown in chapter 11, the essential part of this pro- 
cess is the stretching along the unstable directions of the flow, and in the crud- 
est approximation the dynamics can be modeled by 1 -dimensional expanding 
maps. This observation motivates introduction of a class of particularly simple 
1 -dimensional systems. 
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Figure 25.3: (a) f(x), the full space sawtooth map 
(25.9), A > 2. (b) f(x), the sawtooth map re- 
stricted to the unit circle (25.12), A = 6. 




Example 25.1 Chains of piecewise linear maps. We start by defining the map f on 
the unit interval as 

Ao-{£ + i-a Hffl • A > 2 ' < 259 > 

and then extending the dynamics to the entire real line, by imposing the translation 
property 

fix + h) = fix) + h h 6 Z . (25. 10) 

As the map is discontinuous at x = 1/2, /(1/2) f's undefined, and the x = 1/2 po/nf 
has fo be excluded from the Markov partition. The map is antisymmetric under the 
x-coordinate flip 

/(«) = -/(-*). (25.11) 

so fhe dynamics will exhibit no mean drift; all odd derivatives of the generating function 
(17.11) with respect to f3, evaluated at/3 = 0, will vanish. 

The map (25.9) is sketched in figure 25.3(a). Initial points sufficiently close to 
either of the fixed points in the initial unit interval remain in the elementary cell for one 
iteration; depending on the slope A, other points jump h cells, either to the right or to 
the left. Repetition of this process generates a random walk for almost every initial 
condition. 

The translational symmetry (25.10) relates the unbounded dynamics on the 
real line to dynamics restricted to the elementary cell - in the example at hand, the unit 
interval curled up into a circle. Associated to f (x) we thus also consider the circle map 

fix) = fix) - [fix)] , x = x-[x]€ [0, 1] (25. 12) 

figure 25.3(b), where [■ ■ ■] stands for the integer part (25.2). For the piecewise linear 
map of figure 25.3 we can evaluate the dynamical zeta function in closed form. Each 
branch has the same value of the slope, and the map can be parameterized by a single 
parameter, for example its critical value a = /(1/2), the absolute maximum on the 
interval [0, 1] related to the slope of the map by a - A/2. The larger A is, the stronger 
is the stretching action of the map. 
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As noted in sect. 25.1.1, the elementary cell cycles correspond to either stand- 
ing or running orbits for the map on the full line: we shall refer to h p e Z as the 
jumping number of the p cycle, and take as the cycle weight 



zWlA„|. 



(25.13) 



The diffusion constant formula (25.8) for 1-dimensional maps is 



D = - 



iH 



2 (n) ( ' 



(25.14) 



where the "mean cycle time" is given by (20.25) 



(n) ( = z 



d 1 



dzmz) 



z=l 



Pk 



I a pi ••• A wi 



(25.15) 



and the "mean cycle displacement squared" by (20.28) 



n- 



i 



f dp m i) 



= -T!(- v > 



k (fi Pl + 



+ n Pk ? 



|A P1 ---A 



Pk 1 



(25.16) 



the primed sum indicating all distinct non-repeating combinations of prime cy- 
cles. The evaluation of these formulas for the simple system of example 25.1 will 
require nothing more than pencil and paper. 



Example 25.2 Unrestricted symbolic dynamics. Whenever A is an integer num- 
ber, the symbolic dynamics is exceedingly simple. For example, for the case A = 6 illus- 
trated in figure 25.3(b), the elementary cell map consists of 6 full branches, with uniform 
stretching factor A - 6. The branches have different jumping numbers: for branches 1 
and 2 we have ft = 0, for branch 3 we have h = +1, for branch 4 h = -1, and finally for 
branches 5 and 6 we have respectively h = +2 andh = -2. The same structure reap- 
pears whenever A is an even integer A = 2a: all branches are mapped onto the whole 
unit interval and we have two ft = branches, one branch for which ft = +1 and one for 
which ft = -I, and so on, up to the maximal jump \ft\ = a - 1. The symbolic dynamics 

is thus full, unrestricted shift in 2a symbols {0+, 1+, . . . , (a - 1)+, (a - 1)_ 1-, 0-}, 

where the symbol indicates both the length and the direction of the corresponding jump. 

For the piecewise linear maps with uniform stretching the weight associated 
with a given symbol sequence is a product of weights for individual steps, t sq = t s t q . For 
the map of figure 25.3 there are 6 distinct weights (25.13): 

t\ = t 2 = z/A 

h = e^z/A, U = e- fi z/A, f 5 =e 2/3 z/A, f 6 = e _2 ^z/A. 
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The piecewise linearity and the simple symbolic dynamics lead to the full cancelation 
of all curvature corrections in (20. 7). The exact dynamical zeta function (15.15) is given 
by the fixed point contributions: 



l/?08, z) = 1 - r 0+ - f _ - ■ 



' f (fl-l) + _ f (a-l)- 



a 



1 +^cosh03;) 



(25.17) 



V 7=1 

The leading (and only) eigenvalue of the evolution operator (25.5) is 



s(J3) = logJi 
a 



1 + J] coshOS;) 



A = 2a, a integer. 



(25.18) 



The flow conservation (20. 1 7) sum rule is manifestly satisfied, so s(0) = 0. The first 
derivative s(0)' vanishes as well by the left/right symmetry of the dynamics, implying 
vanishing mean drift (Jc) = 0. The second derivative s(J3)" yields the diffusion constant 
(25.14): 



1 / 9\ 2 l 2 2 2 

n)[ = 2a— = 1, (Jr) = 2— +2— +2— + ■ ■ ■ + 2 
H A \ IC A A A 



)2 ,2 2 2 

'A - V/ f "A +2 A +2 A 
Using the identity 2Li ^ 2 = + 1)( 2 « + l)/6 we obfa/n 



- D 2 
A 



D = ^(A-lXA-2), 



A even integer. 



Similar calculation for odd integer A = 2fc - 1 y/'e/ds 

1 ^ 

£> = — (A 2 - 1) , A odd integer. 



(25.19) 

(25.20) 

exercise 25.1 
(25.21) 



25.2.1 Higher order transport coefficients 

The same approach yields higher order transport coefficients 

S 2 = D, (25.22) 

{3=0 

known for k > 2 as the Burnett coefficients. The behavior of the higher or- 
der coefficients yields information on the relaxation to the asymptotic distribution 
function generated by the diffusive process. Here x t is the relevant dynamical 
variable and S^'s are related to moments (xlfj of arbitrary order. 

Were the diffusive process purely Gaussian 

e tm = | djt ^ e -^D t ) = ^ Dt (2523) 

y4nDt J-oo 

the only St coefficient different from zero would be S2 = D. Hence, nonvan- 
ishing higher order coefficients signal deviations of deterministic diffusion from a 
Gaussian stochastic process. 
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Figure 25.4: (a) A partition of the unit interval 
into six intervals, labeled by the jumping number 
h(x) I = (0 + , 1+, 2+, 2_, 1_, 0_). The partition is 
Markov, as the critical point is mapped onto the 
right border of M\ + . (b) The transition graph for 
this partition, (c) The transition graph in the com- 
pact notation of (25.26) (introduced by Vadim Mo- 
roz). 




Example 25.3 S4 Burnett coefficient. For the map under consideration the first 
Burnett coefficient coefficient £ 4 (or kurtosis (17.21)) is easily evaluated. For example, 
using (25.18) in the case of even integer slope A = 2a we obtain exercise 25.2 



1 



41-60 



(a - l)(2a - l)(4a 2 - 9a + 7) . 



(25.24) 



We see that deterministic diffusion is not a Gaussian stochastic process. Higher 
order even coefficients may be calculated along the same lines. 



25.2.2 Finite Markov partitions 

For piecewise-linear maps exact results may be obtained whenever the critical 
points are mapped in finite numbers of iterations onto partition boundary points, 
or onto unstable periodic orbits. We will work out here an example for which 
this occurs in two iterations, leaving other cases as exercises. The key idea 
is to construct a Markov partition (11.2), with intervals mapped onto unions of 
intervals. 



Example 25.4 A finite Markov partition. As an example we determine a value 
of the parameter 4 < A < 6 for which /(/(1/2)) = 0. As in the integer A case, 
we partition the unit interval into six intervals, labeled by the jumping number h(x) e 
{A1o + , Ati + , Mi+, Mi_, Mi_, Mo_}, ordered by their placement along the unit interval, 
figure 25.4(a). 

In general the critical value a = /(1/2) will not correspond to an interval border, 
but now we choose a such that the critical point is mapped onto the right border of 
Mi + . Equating /(1/2) with the right border of M\ + , x = 1/A, we obtain a quadratic 
equation with the expanding solution A = 2( V2+ 1). For this parameter value f(M\J = 
Mo + {JM\ + , f{Mi_) = Mo_ \JMi_, while the remaining intervals map onto the whole 
unit interval M. The transition matrix (14.1) is given by 



4>' = T<f> 



^ 1 1 1 1 
1110 1 



10 1 
10 1 
10 111 

10111; 



(<f>0 + \ 

4>2- 

01- 



(25.25) 
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One could diagonalize (25.25) on a computer, but, as we saw in chapter 14, the tran- 
sition graph of figure 25.4 (b) corresponding to map figure 25.4 (a) offers more insight 
into the dynamics. Figure 25.4 (b) can be redrawn more compactly as transition graph 
figure 25.4 (c) by replacing parallel lines in a graph by their sum 

— «>>> . = h + h + h . (25.26) 

The dynamics is unrestricted in the alphabet 

J{ = {0 + , 1 + , 2 + + , 2+1+, 2_1_, 2_0_, 1_, 0_) . 

Applying the loop expansion (15.15) of sect. 15.3, we are led to the dynamical zeta 
function 




l/£(J3,z) = 1 - f 0+ - t u - r 2+ o + - t 2+u - t 2 _i_ - t 2 _ Q _ - h_ - t _ 

= 1-^(1+ cosh(/?)) - ^ (cosh(2/3) + cosh(3y6)) . (25.27) 

For grammar as simple as this one, the dynamical zeta function is the sum over fixed 
points of the unrestricted alphabet. As the first check of this expression for the dynam- 
ical zeta function we verify that 



as required by the flow conservation (20. 1 7). Conversely, we could have started by 
picking the desired Markov partition, writing down the corresponding dynamical zeta 
function, and then fixing A by the 1 /£(0, 1) = condition. For more complicated transi- 
tion graphs this approach, together with the factorization (25.35), is helpful in reducing 
the order of the polynomial condition that fixes A. 

The diffusion constant follows from (25. 14) exercise 25.3 

1 2 , 9X l 2 2 2 3 2 
Ml - 4 X +4 A*« H = 2 X +2 A^ +2 A^ 

D = (25.28) 
16 + 8 V2 



It is by now clear how to build an infinite hierarchy of finite Markov partitions: 
tune the slope in such a way that the critical value f(l/2) is mapped into the fixed 
point at the origin in a finite number of iterations p f p ( 1 /2) = 0. By taking higher 
and higher values of p one constructs a dense set of Markov parameter values, 
organized into a hierarchy that resembles the way in which rationals are densely 
embedded in the unit interval. For example, each of the 6 primary intervals can 
be subdivided into 6 intervals obtained by the 2-nd iterate of the map, and for the 
critical point mapping into any of those in 2 steps the grammar (and the corre- 
sponding cycle expansion) is finite. So, if we can prove continuity of D = D(A), 
we can apply the periodic orbit theory to the sawtooth map (25.9) for a random 
"generic" value of the parameter A, for example A = 4.5. The idea is to bracket 
this value of A by a sequence of nearby Markov values, compute the exact diffu- 
sion constant for each such Markov partition, and study their convergence toward 
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the value of D for A = 4.5. Judging how difficult such problem is already for a 
tent map (see sect. 15.5), this is not likely to take only a week of work. 

Expressions like (25.20) may lead to an expectation that the diffusion coeffi- 
cient (and thus transport properties) are smooth functions of parameters control- 
ling the chaoticity of the system. For example, one might expect that the diffusion 
coefficient increases smoothly and monotonically as the slope A of the map (25.9) 
is increased, or, perhaps more physically, that the diffusion coefficient is a smooth 
function of the Lyapunov exponent A. This turns out not to be true: D as a func- 
tion of A is a fractal, nowhere differentiable curve illustrated in figure 25.5. The 
dependence of D on the map parameter A is rather unexpected - even though for 
larger A more points are mapped outside the unit cell in one iteration, the diffusion 
constant does not necessarily grow. 
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Figure 25.6: (a) A map with marginal fixed point, 
(b) The map restricted to the unit circle. 



(a) 




This is a consequence of the lack of structural stability, even of purely hyper- 
bolic systems such as the Lozi map and the 1 -dimensional diffusion map (25.9). 
The trouble arises due to non-smooth dependence of the topological entropy on 
system parameters - any parameter change, no mater how small, leads to creation 
and destruction of infinitely many periodic orbits. As far as diffusion is concerned 
this means that even though local expansion rate is a smooth function of A, the 
number of ways in which the trajectory can re-enter the initial cell is an irregular 
function of A. 

The lesson is that lack of structural stability implies lack of spectral stability, 
and no global observable is expected to depend smoothly on the system parame- 
ters. If you want to master the material, working through one of the deterministic 
diffusion projects on ChaosBook . org/pages is strongly recommended. 



25.3 Marginal stability and anomalous diffusion 

What effect does the intermittency of chapter 24 have on transport properties? A 
marginal fixed point affects the balance between the running and standing orbits, 
thus generating a mechanism that may result in anomalous diffusion. 

Example 25.5 Anomalous diffusion. Consider a 1 -dimensional map of the real line 
on itself shown in figure 25.6 (a), with the same properties as in sect. 25.2, except for a 
marginal fixed point at x = 0. The corresponding circle map is given in figure 25.6 (b). 
As in sect. 24.2.1, a branch with support in Mi, i = 1,2,3,4 has constant slope A,-, 
while f\ Mo is of intermittent form. To keep you nimble, this time we take a slightly 
different choice of slopes. The toy example of sect. 24.2. 1 was cooked up so that the 
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1 Is branch cut in dynamical zeta function was the whole answer. Here we shall take a 
slightly different route, and pick piecewise constant slopes such that the dynamical zeta 
function for intermittent system can be expressed in terms of the Jonquiere funcffemark 25.7 



(25.29) 



k=\ 



Once the fixed point is pruned away the symbolic dynamics is given by the 
infinite alphabet {1, 2,3,4, O'l , 7 2, 0*3, 0'4), i,j,k,l= 1,2,... (compare with table 24.1). 
The partitioning of the subinterval Mo is induced by M&( r i g ht) = fhtht)^ 3 U-^U) 
(where f^ gM) denotes the inverse of the right branch of /\m ) and the same reason- 
ing applies to the leftmost branch. These are regions over which the slope of /U< is 
constant. Thus we have the following stabilities and jumping numbers associated to 
letters: 



0*3, 0*4 




q/2 


h p - 1 




O'l, 0'2 


Ap 


- q/2 


h p - -1 




3, 4 


A, 


= +A 


hp = 1 




2, 1 


A P 


= +A 


hp = -1 , 


(25.30) 



where a = l/s is determined by the intermittency exponent (24.1), while q is to be 
determined by the flow conservation (20. 1 7) for f: 



A 

(where ( is the Riemann zeta function), so thatq = (A-4)/(2A^(a+ 1)). The dynamical 
zeta function picks up contributions just by the alphabet's letters, as we have imposed 
piecewise linearity, and can be expressed in terms of a Jonquiere function (25.29): 



l/io(z,P) = l--zcosh^- -A 4 - zcosh/3-J(z,a + l). 
A Af (1 + a) 

Its first zero z(J3) is determined by 

4 A-4 1 
—z + : — — -z ■ J(z, a + 1) = 



(25.31) 



A A£(l + or) 



COShyS 



D vanishes by the implicit function theorem, z"(fi)\t 



when a < 1. The 



physical interpretation is that a typical orbit will stick for long times near the 
marginal fixed point, and the 'trapping time' will be larger for higher values of 
the intermittency parameter s (recall a = l/s). As always, we need to look more 
closely at the behavior of traces of high powers of the transfer operator. 

The evaluation of transport coefficient requires one more derivative with re- 
spect to expectation values of state space observables (see sect. 25.1): if we use 
the diffusion dynamical zeta function (25.7), we may write the diffusion coeffi- 
cient as an inverse Laplace transform, in such a way that the distinction between 
maps and flows has vanished. In the case of 1 -dimensional diffusion we thus have 



D = lim 

t—>oo 



d 2 i r 

d/3 2 2ni J a . 



° st C(fi,s) 
ase 

m s) 



(25.32) 



/?=o 
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where the £' refers to the derivative with respect to s. 

The evaluation of inverse Laplace transforms for high values of the argument 
is most conveniently performed using Tauberian theorems. We shall take 



-r 

Jo 



co(A) =| dxe Ax u(x), 



with u(x) monotone as x —> oo; then, as A i-» and x i-» oo respectively (and 
P e (O.oo), 

if and only if 

w(x) ~ xP~ l L(x) , 

T(p) 

where L denotes any slowly varying function with lim,_>oo L(ty)/L(t) = 1. Now 

l/Zo(e-',fi) 1 - cosh ,6 - x^ye"*(e" s , a + 1) coshpV ' 
Taking the second derivative with respect to /? we obtain 
J 2 



d02 (i/^v^/rV^ 



^ + Jnh(J(e- s ,a + l) + J(e- s ,a)) 

g a (s). (25.33) 



(i-x^-a^-^>«+D) 



The asymptotic behavior of the inverse Laplace transform (25.32) may then be 
evaluated via Tauberian theorems, once we use our estimate for the behavior of 
Jonquiere functions near z = 1 - The deviations from normal behavior correspond 
to an explicit dependence of D on time. Omitting prefactors (which can be calcu- 
lated by the same procedure) we have 

( s~ 2 for a > 1 

ga(s) ~ s~ (a+1) for a e (0, 1) 
{ l/(s 2 lns) for ar= 1. 

The anomalous diffusion exponents follow: exercise 25.6 

P for a e (0, 1) 
({x-x ) 2 ) t ~ \ t/\nt for a= 1 (25.34) 
t for a > 1 . 
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Resume 

Perfection itself is imperfection. 
— Vladimir Horowitz 

With initial data accuracy 8x = |5x(0)| and system size L, a trajectory is predictable 
only to the finite Lyapunov time 7Ly ap * A ln\L/6x\ . Beyond that, chaos rules. 
We have discussed the implications in sect. 1.8: chaos is good news for prediction 
of long term observables such as transport in statistical mechanics. 

The classical Boltzmann equation for evolution of 1 -particle density is based 
on stosszahlansatz, neglect of particle correlations prior to, or after a 2-particle 
collision. It is a very good approximate description of dilute gas dynamics, but 
a difficult starting point for inclusion of systematic corrections. In the theory 
developed here, no correlations are neglected - they are all included in the cycle 
averaging formula such as the cycle expansion for the diffusion constant 

n = 1 1 Y'r + "- + / W 2 

2d(T)^ K lA^—A^I ■ 

Such formulas are exact; the issue in their applications is what are the most ef- 
fective schemes of estimating the infinite cycle sums required for their evaluation. 
Unlike most statistical mechanics, here there are no phenomenological macro- 
scopic parameters; quantities such as transport coefficients are calculable to any 
desired accuracy from the microscopic dynamics. 

For systems of a few degrees of freedom these results are on rigorous footing, 
but there are indications that they capture the essential dynamics of systems of 
many degrees of freedom as well. 

Though superficially indistinguishable from the probabilistic random walk 
diffusion, deterministic diffusion is quite recognizable, at least in low dimen- 
sional settings, through fractal dependence of the diffusion constant on the system 
parameters, and through non-Gaussion relaxation to equilibrium (non-vanishing 
Burnett coefficients). 

That Smale's "structural stability" conjecture turned out to be wrong is not a 
bane of chaotic dynamics - it is actually a virtue, perhaps the most dramatic ex- 
perimentally measurable prediction of chaotic dynamics. As long as microscopic 
periodicity is exact, the prediction is counterintuitive for a physicist - transport 
coefficients are not smooth functions of system parameters, rather they are non- 
monotonic, nowhere differentiable functions. 

Actual evaluation of transport coefficients is a test of the techniques developed 
above in physical settings. In cases of severe pruning the trace formulas and er- 
godic sampling of dominant cycles might be more effective strategy than the cycle 
expansions of dynamical zeta functions and systematic enumeration of all cycles. 
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Commentary 

Remark 25.1 Lorentz gas. The original pinball model proposed by Lorentz [25.4] 
consisted of randomly, rather than regularly placed scatterers. 

Remark 25.2 Who's dunnit? Cycle expansions for the diffusion constant of a particle 
moving in a periodic array have been introduced by R. Artuso [25.5] (exact dynamic- 
al zeta function for 1-dimensional chains of maps (25.8)), by W.N. Vance [25. 6], and by 
P. Cvitanovic, J. -P. Eckmann, and P. Gaspard [25.7] (the dynamical zeta function cycle 
expansion (25.8) applied to the Lorentz gas). 

Remark 25.3 Lack of structural stability for D. Expressions like (25.20) may lead 
to an expectation that the diffusion coefficient (and thus transport properties) are smooth 
functions of the chaoticity of the system (parameterized, for example, by the Lyapunov 
exponent A = In A). This turns out not to be true: D as a function of A is a fractal, 
nowhere differentiable curve shown in figure 25.5. The dependence of D on the map 
parameter A is rather unexpected - even though for larger A more points are mapped 
outside the unit cell in one iteration, the diffusion constant does not necessarily grow. 
We refer the reader to refs. [25.15, 25.16] for early work on the deterministic diffusion 
induced by 1-dimensional maps. The sawtooth map (25.9) was introduced by Grossmann 
and Fujisaka [25.17] who derived the integer slope formulas (25.20) for the diffusion 
constant. The sawtooth map is also discussed in refs. [25.18]. The fractal dependence of 
diffusion constant on the map parameter is discussed in refs. [25.9, 25.8, 25.10]. Sect. 1.8 
gives a brief summary of the experimental implications; for the the current state of the art 
of fractal transport coefficients consult the first part of Klage's monograph [25. 12]. Would 
be nice if someone would eventually check these predictions in experiments... Statistical 
mechanicians tend to believe that such complicated behavior is not to be expected in 
systems with very many degrees of freedom, as the addition to a large integer dimension 
of a number smaller than 1 should be as unnoticeable as a microscopic perturbation of a 
macroscopic quantity. No fractal-like behavior of the conductivity for the Lorentz gas has 
been detected so far [25.14]. (P. Cvitanovic and L. Rondoni) 

Remark 25.4 Symmetry factorization in one dimension. In the p = limit the 
dynamics (25.1 1) is symmetric under x — » -x, and the zeta functions factorize into prod- 
ucts of zeta functions for the symmetric and antisymmetric subspaces, as described in 
sect. 21.1.1: 



The leading (material flow conserving) eigenvalue z — 1 belongs to the symmetric sub- 
space l/f s (0, 1) = 0, so the derivatives (25.15) also depend only on the symmetric sub- 



1 



1 



1 



dz( 



1 <3 1 15 1 

f s 8z £„ £j dz £s ' 



(25.35) 



space: 



{n)r = 



d 1 



7 



~dzt(0,z) c 



= i 
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length 


# cycles 


«0,0) 


A 


1 


5 


-1.216975 




2 


10 


-0.024823 


1.745407 


3 


32 


-0.021694 


1.719617 


4 


104 


0.000329 


1.743494 


5 


351 


0.002527 


1.760581 


6 


1243 


0.000034 


1.756546 



Table 25.1: Fundamental domain, w=0.3 . 



_J d_ 1_ 

U^z) Z dl £(0,z) t _ 

Implementing the symmetry factorization is convenient, but not essential, at this level of 
computation. 

Remark 25.5 Lorentz gas in the fundamental domain. The vector valued nature 
of the generating function (25.3) in the case under consideration makes it difficult to 
perform a calculation of the diffusion constant within the fundamental domain. Yet we 
point out that, at least as regards scalar quantities, the full reduction to M. leads to better 
estimates. A proper symbolic dynamics in the fundamental domain has been introduced 
inref. [25.19]. 

In order to perform the full reduction for diffusion one should express the dynamical 
zeta function (25.7) in terms of the prime cycles of the fundamental domain M. of the 
lattice (see figure 25.2) rather than those of the elementary (Wigner-Seitz) cell M. This 
problem is complicated by the breaking of the rotational symmetry by the auxiliary vector 
P, or, in other words, the non-commutativity of translations and rotations: see ref. [25.7]. 

Remark 25.6 Anomalous diffusion. Anomalous diffusion for 1-dimensional inter- 
mittent maps was studied in the continuous time random walk approach in refs. [24.10, 
24. 1 1]. The first approach within the framework of cycle expansions (based on truncated 
dynamical zeta functions) was proposed in ref. [24.12]. Our treatment follows methods 
introduced inref. [24.13], applied there to investigate the behavior of the Lorentz gas with 
unbounded horizon. 



Remark 25.7 Jonquiere functions. In statistical mechanics Jonquiere function 
(25.29) appears in the theory of free Bose-Einstein gas, see refs. [24.22, 24.23]. 



(25.36) 
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Exercises 



25.1. Diffusion for odd integer A. Show that when the 
slope A = 2k - 1 in (25.9) is an odd integer, the dif- 
fusion constant is given by D = (A 2 - l)/24, as stated in 
(25.21). 

25.2. Fourth-order transport coefficient. Verify (25.24). 
You will need the identity 

" 1 

^ k 4 = —n(n + 1)(2« + l)(3n 2 + 3« - 1) . 
k=\ 



versa, so that 1, 2, 3, 4 are standing while leads to both 
positive and negative jumps. Build the corresponding 
dynamical zeta function and show that 



for a > 2 
for a - 2 
f- u for a e (1,2) 
f 2 /lnf for a — 1 
f- for a e (0, 1) 



/ 

tint 

3-a 



25.3. Finite Markov partitions. 



25.4. Maps with variable peak shape: 

lowing piecewise linear map 



25.5. 



Verify (25.28). 

Consider the fol- 



3x 

is 



Mx) = 



!~(§l^-*l) 
i-M x -k 2+s >) 



x e Mi 

xeM 2 (25.37) 
x e M3 



where Mi = [0, ±(1 - 6)], M 2 = [f(l - 5), 1(2 + 5)], 

M3 = [g(2 + 6), ij, and the map in [1/2, 1] is obtained 
by antisymmetry with respect to x — 1 /2, y = 1 /2, Write 
the corresponding dynamical zeta function relevant to 
diffusion and then show that 



D = 



6(2 + 6) 
4(1-5) 



See refs. [25.21, 25.22] for further details. 

Two-symbol cycles for the Lorentz gas. Write down 
all cycles labeled by two symbols, such as (0 6), (1 7), 
(1 5) and (0 5). 

ChaosBook.org/pages offers several project-length 
deterministic diffusion exercises. 



25.6. Accelerated diffusion, (medium difficulty) Consider 
a map h, such that h — f of figure 25.6 (b), but now run- 
ning branches are turned into standing branches and vice 



25.7. Recurrence times for Lorentz gas with infinite hori- 
zon. Consider the Lorentz gas with unbounded 
horizon with a square lattice geometry, with disk ra- 
dius R and unit lattice spacing. Label disks accord- 
ing to the (integer) coordinates of their center: the se- 
quence of recurrence times \tj] is given by the set of 
collision times. Consider orbits that leave the disk sit- 
ting at the origin and hit a disk far away after a free 
flight (along the horizontal corridor). Initial conditions 
are characterized by coordinates ((/>, a) (0 determines the 
initial position along the disk, while a gives the an- 
gle of the initial velocity with respect to the outward 
normal: the appropriate measure is then dip cos a da 
(<p € [0, 27r), a 6 [-n/2,n/2]. Find how <p(T) scales 
for large values of T: this is equivalent to investigating 
the scaling of portions of the state space that lead to a 
first collision with disk (n, 1), for large values of n (as 
11 h 00 n ^ T). 

25.8. Diffusion reduced to the fundamental domain. 



Maps such as figure 25.3 are antisymmetric. Re- 
duce such antisymmetric maps as in example 9.4, and 
write down the formula (25.14) for the diffusion con- 
stant D in terms of the fundamental domain cycles (rela- 
tive periodic orbits) alone (P. Gaspard says it cannot be 
done [25.7]). 
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